Moment Analysis of Multiplicity 
Distributions 



A. Capella, I.M. DreminQ, V.A. Nechitailo\ J. Tran Thanh Van 

Laboratoire de Physique Theorique et Hautes Energies, 
Bat. 211, Universite Paris-Sud, 91405 Orsay, France 

Abstract 

Moment analysis of global multiplicity distribution previously done 
for e'^e~ and hh processes is applied to hA and AA collisions. The 
oscillations of cumulants as functions of their rank are found in all 
the cases. Some phenomenological approaches and quark-gluon string 
models are confronted to experimental data. It has been shown that 
the analysis is a powerful tool for revealing the tiny features of the 
distributions, and its qualitative results in various processes are rather 
stable for different multiplicity cut-offs determined by experimental 
(or Monte-Carlo) statistics. 

1 Introduction 

Multiplicity distributions of multiparticle reactions contain in the integrated 
form all the correlations of the system. At the same time they are measured 
with best accuracy in experiment. Therefore their study is of particular 
interest. The shapes of the distributions differ drastically for various pro- 
cesses at different energies. Phenomenologically, it is popular to use Poisson 
(or even sub-Poisson) distributions at lower energies and negative binomial 
distribution (nowadays, modified negative binomial distribution) at higher 
energies. The parameters of the distributions vary for different cases and are 
not determined with high enough precision when fits are done. 

However, any multiplicity distribution can be represented not only by 
the probabilities of n-particle events but also by its moments or by its gen- 
erating function. It happens that such description is preferred both from 
the theoretical point of view and for the stability of results obtained from 
experimental data for various reactions. The solution of QCD equations for 
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generating functions predicts |1|, 0] the oscillatory behaviour of the cumulant 
moments of the parton distribution in QCD jets as functions of their rank 
with first minimum located at the point determined by the inverse power of 
the anomalous dimension 7. 

Surprisingly enough, the oscillations of the ratio of cumulant to factorial 
moments with similar periodicity have been found in experimental hadron 
multiplicity distributions in e+e~ and pp{pp) collisions at high energies 
No pattern of this kind appears in phenomenological fits (for review see 0). 

In this paper we apply the moment analysis to multiplicity distributions 
in hadron-nucleus and nucleus-nucleus collisions. The similar (to ee and 
hh) shape of the above ratio appears in both cases. The only difference is 
in the depth of the first minimum for different projectiles and targets but 
not in its location. We compare these results to purely phenomenological 
fits of the distributions and to some theoretical models (DPM - dual parton 
model and QGSM - quark-gluon string model). While phenomenological fits 
do not describe shapes of the ratio obtained from experimental data, the 
string models reproduce them qualitatively rather well. The minor quanti- 
tative difference can be accounted for by slight variation of parameters and 
should be attributed to experimental uncertainties that demonstrate also how 
sensitive the method of moment analysis is to tiny features of multiplicity 
distributions. In that aspect, the qualitative stability of the ratio shapes 
for reactions with drastically different distributions looks very impressive. 
The general feature of branching processes inherent both in QCD and in 
theoretical Monte-Carlo models can be in charge of it. 

We have checked that the truncation of multiplicity distributions due 
to finite statistics and conservation laws does not influence drastically this 
dependence even though it produces some additional oscillation effect which 
should be taken into account when comparing models to experimental data. 
The theoretical arguments in favor of it can be got from studies in paper [^], 
where it has been shown that conservation laws give rise to corrections of 
the order of 7^, while the oscillation effect appears already due to modified 
leading logarithm approximation terms of the order of 7^. 

Besides, we have found the zeros of the truncated generating functions. 
Similar to ee and hh processes, they tend to lie close to the circle of unit 
radius in the complex plane for the high enough truncation limit. Since the 
highest multiplicities are in AA-collisions, the circle is most clearly seen there. 
The rightmost zeros tend to the real axis at higher multiplicities that reminds 
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of Lee- Yang conjecture in statistical physics about the phase transition point 
to which the rightmost zeros of the grand partition function approach when 
volume increased. It implies that the singularity of the generating function 
is very close to the point where all the moments are calculated as derivatives 
of the generating function. That explains why the moment analysis is so 
sensitive to slight variations of multiplicity distributions and, at the same 
time, provides qualitatively similar results for various processes. 



2 Moment Analysis of Hadron-Nucleus and 
Nucleus-Nucleus Collisions 

First, let us describe the generalities of the moment analysis. The normalized 
factorial (Fg) and cumulant {Kg) moments of the multiplicity distribution P„ 
are defined as 



(n(n-l)...(n-g + l)) _ ^ n(n - 1) . . . (n - g + 1)P„ 
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where P„ is a probability of ra-particle events, q is the rank of the moment, 
= n\ / m\{n — m)\ are the binomial coefficients. 

Thus, if the multiplicity distribution P„ is known one calculates the fac- 
torial moments of any rank according to and, afterwards, using the 
recurrence relations finds out the cumulants. 

Theoretically, it is more convenient to start with the generating function 



n=0 



and calculate the moments as its derivatives 

1 d^G 
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In practice, the multiplicity distribution is known up to some maximum 
multiplicity N, and one has to deal with the truncated generating function 

N 

Gn = J2 '^"^n. (6) 
n=0 

The interest to the moment analysis of multiplicity data appeared after the 
solution of QCD equations for the generating function was found 0]- It 
predicts a special oscillation pattern for the ratio of cumulant to factorial 
moments 

H. = ^ (7) 

with first minimum located at qmin ~ 5 and determined by the inverse value 
of the QCD anomalous dimension 

70 = {2N,as/7cy/\ (8) 

where as is a coupling constant, A^^c = 3 is the number of colours. For further 
details see papers || and the review paper 0. 

Surely, it was the prediction for the moments of parton (mostly, gluon) 
multiplicity distributions. However, when applied to final hadrons in e^e~ 
and pp{pp) experimental events the analysis shows the similar structure 
of the ratio Hq. 

We extend this analysis to hadron-nucleus and nucleus-nucleus collisions 
using experimental data, theoretical Monte-Carlo models and phenomeno- 
logical fits. We check also how important are various multiplicity cut-offs for 
the moments. The zeros of the truncated generating function at various 
cut-offs are found and discussed. 

This analysis is of a special interest because the shapes of multiplicity 
distributions themselves differ strongly from those in ee and pp collisions. 

In Fig.l we plot the ratio Hg calculated according to the dual parton 
model and quark-gluon string model for hadron-nucleus and nucleus- 
nucleus collisions. For a comparison, we show this ratio for charged particles 
in proton-antiproton collisions at 546 GeV calculated from experimental data 
and published in [Q. 

One concludes that the general shapes of Hg are similar for all these 
processes. The oscillations with first negative minimum in the range g = 4— 6 
are clearly seen. Its depth increases at each step from hh to hA and AA. 
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The stability of qualitative features of moments is especially remarkable if 
one compares the multiplicity distributions in all the cases which are very 
different, indeed. 

No phenomenological fit, we are aware of, is able to reproduce such a 
behaviour. For Poisson distribution all cumulants (and, consequently, ratios 
Hq) are identically equal to zero (at q > 2). The modified negative binomial 
distribution has the generating function 

with three free parameters r. A, k related to the average multiplicity by (n) = 
k{r — A). The negative binomial distribution is obtained for A = 0. Using 
(H) it is easy to show that 

Since k > 0, the cumulants are always positive for the negative binomial 
distribution. Recent fits of e'^e" data H, p| give rise to A < and r < 



|A| < 1. According to (|l^) it would imply that cumulants change sign at 
each q being negative at even values of q and positive at odd ones. They do 
not reveal the broad oscillations shown in Fig.l. The whole picture reminds 
somewhat the pattern of fixed multiplicity distribution (see 0]). In general, 
we would like to comment that the fits used in look rather strange since 
they imply sub-Poissonian distribution of negatively charged particles (with 
K2 < 0), while it is usually claimed that it is super-Poissonian one at 
high energies. Moreover, with the set of parameters A = — 0.75,r = 0.65 
available in [§, ^ one gets negative probabilities for large n as is easily seen 
from d^) and was confirmed by computer calculations. Moment analysis 
of charged particle distributions in e+e" data done in |10] has shown 



the pattern reminding that of pp and, according to Fig.l, of reactions with 
nuclei. The only qualitative difference is the smaller depth of the minimum 
in ee compared to hadronic reactions. 

Contrary to phenomenological fits, the dual parton model and the 
quark-gluon string model look quite successful in describing the qualitative 
behaviour of the ratio Hg in hadron-nucleus collisions as shown in Figs. 2-4, 
where the experimental data of NA22 collaboration at 250 GeV have been 
used W^. The general periodicity of the curves and the absolute values at 
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minima and maxima are reproduced rather well even though their positions 
are somewhat shifted for pAl and KAl. There are several factors which 
influence the shapes of iJ^-curves, and we discuss them in the next section. 



3 Asymptotically Subleading Terms of the Mo- 
ments 



The stability of the oscillating pattern in Fig.l is astonishing because, first, 
the pattern differs from those of widely used in probability theory and, sec- 
ond, the multiplicity distributions in various processes under consideration 
are very different and do not seem to have much in common (contrary to 
ifg-ratios). Moreover, this pattern seems sometimes very sensitive to tiny 
modifications of multiplicities themselves. This is related to the subtraction 
procedure used to calculate cumulants according to eq. @ when factorial 
moments are known. 

However, when comparing to experimental data one should be sure that 
the proper distributions have been chosen. By that, for example, we mean 
that the moment analysis of charged multiplicities provides different values 
of ratios H^'^'^^ compared to the values -f^^^-* obtained for negatively charged 
particles. If one considers e"^e~ collisions, the number of charged particles 
is twice that of negatives. Therefore, the corresponding generating functions 
are related by 




(11) 



Herefrom, one gets for the first five moments: 




(12) 
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and the same relations are valid for cumulants (one should just replace F by 
in (|12D). One concludes that F^"^^ > F^'^ at any g, while the inequality 

xlf^^ > K2 ^ is always valid for cumulants but K^^'^^ can be less than ^ 
if the distribution of negatives is sub-Poiisonian one (i.e., if < 0). If 
negative particles are distributed according to Poisson law i.e. Kj^"^ = for 
g > 2, the charged particles dispersion differs from Poissonian. In that case, 
K2^^'^ = {n)~^ > but higher cumulants are equal to zero according to (|1^) 
as if Poisson law is restored just for them. 

One can think about this procedure as if the negative particles are created 
by a "cluster" of negative+positive particle (due to charge conservation they 
are always produced in pairs). In general, one can consider a cluster model 
with a definite probability for cluster production and its decay into k particles 
(see for example). In that case the relation between the generating 
functions for particles and clusters is 

Gp{z) = Ge(^'), (13) 

and the formulae (|T2[) can be easily generalized to 



k 



piP) _ p(c) I 3(^-l) p(c) I {k-l){k-2) 

etc. Herefrom, one gets particle distributions wider than cluster ones because 
^2^"* > ^2^'' for k > 1, e.g., super-Poissonian ones for particles created by 
Poisson clusters typical for asymptotical multiperipheral and string models. 
In particular, for the model of Poissonian clusters with k = 1.4 charged 
particles per cluster 0] one gets 



X?) = > 0; = < 0. (15) 

To get K^^^ > 0, one needs k > 2. 

If KNO-scaling holds at asymptotically high energies, all Fg tend to con- 
stants and (rich) 00. Therefore, all the moments of both distributions 
coincide at asymptotics. At present energies, however, the correction terms 
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of the order 0(1/ (nch)) can be still important, especially, for those cumu- 
lants which are close to zero. Therefore, cumulant analysis can give different 
results for different distributions. Above, we apphed it to charged particles 
multiplicities. 

The finite energy of colliding particles and the final experimental statis- 
tics truncate the multiplicity distribution at some finite multiplicity A^. The 



relations (p!2D, (|Tj) are independent of this truncation since they are valid for 
truncated generating functions (||) as well as for total ones (H). However, it 
influences the values of the moments in these formulae. E.g., the truncated 
Poisson law for negative particles would give rise to non-zero values of cu- 



mulants. Actually, it was conjectured in ||13| , |14| that at present energies this 
effect could be rather important for ee and pp collisions since it imitates the 
oscillations of cumulants. Surely, it should disappear at higher energies. This 
is clear both from model calculations and from theoretical arguments of 
where this effect is of the order of 7^, while the oscillations in QCD appear 
already at the modified logarithm approximation of 7^ order [|l], ^. Never- 
theless, at present energies it should be carefully treated for each experiment 
together with uncertainties imposed by the error bars due to finite statistics. 
In any Monte-Carlo model they can be reduced by enlarged statistics which 
is a matter of computer time only. 

We do not show the error bars in our Figs but the stability of curves and 
the regularity of trends shows that they hardly can change our conclusions. 
For ee processes the minima are smaller but even there it was shown |lOl 
that the first minimum is still reliable. We show in Fig. 5 that the shapes 
of ifg-curves in nucleus-nucleus collisions are not very sensitive to the mul- 
tiplicity cut-off. Truncating the rather fiat multiplicity distribution at two 
different values of A^ we do not observe any strong shift of the curve and 
noticeable change of its shape. 



4 Zeros of the Truncated Generating Func- 
tion 

The truncated generating function (^) is a polynomial of A^-th power with 
positive coefficients. Therefore, it possesses A^ complex conjugate zeros in 
the complex z-plane. We have found their locations in all the above cases 
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varying the value of N. 

The general picture becomes quite stable at large A^. The zeros tend to 
lie close to unit circle in z-plane. The rightmost ones move to the real axis 
with increasing as if trying to pinch it at some value of z slightly exceeding 
1. They move closer and closer to 1 if we compare ee, hh, hA, AA collisions, 
correspondingly. Their convergence point implies that the singularity of the 
total generating function is located near z = 1, where all the moments are 
calculated (see eqs. (|),®). It means that moments are very sensitive to the 
location and the origin of the singularity. (Let us note that the singularity of 
the negative binomial distribution appears at z = 1 + according to (|^), 
i.e. quite far from 2; = 1 if r < 1 as proposed in that contradicts to Figs 
below) . 

We demonstrate all these features in Figs. 6-9. To show how close are 
the rightmost zeros to the point 2; = 1 in different processes we give their 
coordinates x,y and the radii R = (x^ + y'^Y^'^ in the Table. Also shown 
are the maximum multiplicities A^. All the radii exceed slightly 1, and the 
nearest ones to 1 are in nucleus-nucleus collisions. 

The above findings could be interesting also if one speculates about the 
statistical analogies in particle physics [^] in terms of Feynman- Wilson liquid 



1^, [13] • Then the truncated generating function reminds the grand partition 
function and the variable z plays the role of fugacity. It has been proven by 
Lee and Yang that its zeros should lie on the circle in the complex plane and 
the tendency of rightmost zeros to pinch the positive real axis means that 
there is a phase transition in the system. Let us note, however, that there 
is no symmetry in P„ which helped Lee and Yang to prove their statement 
about zeros location. The similar patterns appearing in all processes in 
particle physics are quite encouraging for such analogies. The nature of the 
singularity could be revealed from analysis of experimental data, in principle. 



5 Discussion and Conclusions 

The moment analysis of multiplicity distributions as applied to hA and AA 
collisions reveals the oscillation pattern of the ratio Hq similar to that previ- 
ously found for ee and hh processes. However, the amplitude of oscillations 
increases for heavier colliding objects, while their periodicity remains stable. 
Phenomenological fits do not show anything similar to observed patterns. 
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The string models reproduce them rather well. The multiplicity cut-off is 
not very important until one comes too close to average multiplicity. The 
zeros of the truncated generating function tend to be near the unit circle in 
the complex plane with rightmost zeros moving at high multiplicities toward 
the real positive axis and pinching it near the point z = 1 where derivatives 
are taken when moments arc calculated. This is the singularity point of the 
total generating function and it is closer to ^ = 1 for heavier colliding objects. 

These findings show that there is much more similarity in different col- 
lision processes than it could be envisaged from rather different multiplicity 
distributions. Probably, it is due to the common branching origin of them. 

Also, the singularity point of the total generating function would be an in- 
teresting topic to learn more about. The density of zeros and their movement 
to the real axis at higher multiplicities in various processes could indicate its 
nature. 

Since the moments are calculated as derivatives of the generating function 
at the point very close to the singularity, the moment analysis is powerful 
and sensitive tool of revealing its structure and tiny features of multiplicity 
distributions. 
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Figure Captions 

Fig.l The ratio Hq for pp, hA and AA collisions. (See text for more details). 

Fig. 2 The ratio Hq for pAl collisions at 250 GeV as calculated from experi- 
mental data |ll| and according to DPM and QGSM 0. 



Fig. 3 The ratio Hq for ttAI collisions at 250 GeV (experimental data [jTT 
and QGSM W). 



Fig. 4 The ratio Hq for i^'Al collisions at 250 GeV (experimental data [jTT 
and QGSM W). 



Fig.5 The ratio Hq for ■^^S^^^U collisions at 200A GeV of •^^S in the labora- 
tory system, calculated from DPM multiplicity distributions truncated 
at two values of maximum multiplicity N. 

Fig. 6 Zeros of the truncated generating function for pAl collisions at 250 
GeV. 

Fig. 7 The same as Fig. 6 but for vrAl. 
Fig. 8 The same as Fig. 6 but for i^Al. 

Fig. 9 Zeros of the truncated generating function for 32g238-[j collisions at 
200A GeV of '^^S in the laboratory system according to DPM 

Table 

Location of the rightmost zeros of the truncated generating functions for 
different processes. 



collision 


R 


y 


X 


N 


vr Al (QGSM) 


1.16042 


0.24190 


1.1349 


41 


TT Al (exp) 


1.16042 


0.24771 


1.1337 


41 


K Al (QGSM) 


1.12553 


0.28954 


1.0877 


32 


K Al (exp) 


1.14096 


0.30445 


1.0996 


32 


p Al (DPM) 


1.15537 


0.29461 


1.1172 


27 


p Al (QGSM) 


1.09168 


0.28592 


1.0536 


27 


p Al (exp) 


1.03598 


0.31132 


0.9881 


27 


S U (DPM) 


1.02265 


0.07061 


1.0202 


80 
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